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ABSTRACT 

We investigate various astrophysical contributions to the statistical uncertainty of precision radial 
velocity measurements of stellar spectra. We first analytically determine the intrinsic uncertainty in 
centroiding isolated spectral lines broadened by Gaussian, Lorentzian, Voigt, and rotational profiles, 
finding that for all cases and assuming weak lines, the uncertainty is the line centroid is ay ~ 

C 0^/^/(IU/q^^), where 0 is the full-width at half-maximum of the line, W is the equivalent width, 
and Jo is the continuum signal-to-noise ratio, with C a constant of order unity that depends on the 
specific line profile. We use this result to motivate approximate analytic expressions to the total radial 
velocity uncertainty for a stellar spectrum with a given photon noise, resolution, wavelength, effective 
temperature, surface gravity, metallicity, macroturbulence, and stellar rotation. We use these relations 
to determine the dominant contributions to the statistical uncertainties in precision radial velocity 
measurements as a function of effective temperature and mass for main-sequence stars. For stars more 
massive than ~ I.l Mq we find that stellar rotation dominates the velocity uncertainties for moderate 
and high resolution spectra {R > 30,000). For less massive stars, a variety of sources contribute 
depending on the spectral resolution and wavelength, with photon noise due to decreasing bolometric 
luminosity generally becoming increasingly important for low-mass stars at fixed exposure time and 
distance. In most cases, resolutions greater than 60,000 provide little benefit in terms of statistical 
precision, although higher resolutions would likely allow for better control of systematic uncertainties. 

We find that the spectra of cooler stars and stars with higher metallicity are intrinsically richer in 
velocity information, as expected. We determine the optimal wavelength range for stars of various 
spectral types, finding that the optimal region depends on the stellar effective temperature, but for 
mid M-dwarfs and earlier the most efficient wavelength region is from 6000 A to 9000 A. 


1. INTRODUCTION 

Current generation radial velocity (RV) surveys for ex¬ 
oplanets mostly focus on relatively bright, V < 8.5, stars 
using single-object spectrographs. Since they are using 
single-object spectrographs, the surveys target one star 
at a time and expose up to a desired continuum signal- 
to-noise ratio (SNR). The V cut-off used by the RV sur¬ 
veys is generally driven by the fain t limit of the Hippar- 
cos results fe.g.. lMarcv et al.ll2005ll : current RV surveys 
deliberately pre-select all of their stars to ensure that 
they will be dwarfs, chromospherically quiet, and do not 
have known stellar companions - so that they will be 
good RV targets. Initially the surveys largely restricted 
their targets to single G- and K-dwarfs ( e.g. the stel¬ 
lar sample described in lWright et al.l[2004ll . This choice 
was driven on the one side by the faintness of M-dwarfs, 
which makes getting a high SNR in a reasonable exposure 
time difficult, and on the other side by the rapid ro ta- 
tion of stars hotter than the Kraft Break (|KraftJII970ll at 
6250K, which widens the stellar lines and makes precision 
RV measurements difficult. Over the last two decades, 
RV surveys have thus surveyed effectively all of the sin¬ 
gle G- and K-dwarfs brighter than V < 8.5 for Jupiter- 
mass planets out to periods of 5.5 years (see, for example, 
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iCumming et al.l[2008t iWrightI[200511 . 

In the last ten years, a new generation of RV sur¬ 
veys arose that t argeted M-dwarf s (M2K and the Kec k 
M-dwarf Survey, IAdds et ^120101 : 1.Tohnson et al. 201C), 
fainter high-metallicity stars fN2K nFischer et~ 2005 1. 
and for mer A-stars that are now sub-giants (“Retired A 
Stars,” [Johnson et aLll^OOTll - to name a few examples. 
These newer surveys use the same observing mode as pre¬ 
vious surveys: single-object spectrographs and exposing 
to a desired SNR. They are nearly complete for Jupiter- 
mass planets out to several hundred days, though they 
have yet to survey all of the possible stars available. 

As the field moves forward, next generation of RV sur¬ 
veys will target fainter stars in broader observing modes. 
This will be enabled by the results from the GAIA mis¬ 
sion, w hich will allow fo r vetting of target stars up to 
V < 20 (|de Bruiinel l20'Tl. and also by the exhaustion of 
unsurveyed bright stars. There will also be a consider¬ 
able demand for RV resources to follow-up planet can¬ 
didates from the upcoming Transit ing Exoplanet Survey 
Satellite ITESS. [Ricker et ffi]l2015i l ar id Planetary Tran- 
sits and Oscillations of stars iPLATO. lRauer et al.ll20I3[l 
mission s, as well as residua .1 planet candidates from the 
Kepler ([Borucki et al.l[2M^ mission. Importantly, these 
photometric exoplanet surveys do not pre-select their 
targets in the same manner as the traditional RV sur¬ 
veys (if they pre-select at all). Vetting the candidates 
from these missions will therefore necessitate precision 
RV measurements of stars that are potentially more ac¬ 
tive, hotter, more evolved, and rotating faster than the 
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targets of the first-generation surveys. 

In all of this, understanding the sources of velocity 
uncertainty in stellar RV measurements is critical. This 
allows us to both appropriately understand the detection 
sensitivities of current RV surveys, and more efficiently 
design and execute future precision RV searches for ex¬ 
oplanets. In particular, consider a multi-object preci¬ 
sion RV survey, which would have design considerations 
very different from present single-star searches. Firstly, 
a multi-object survey will image a much narrower spec¬ 
tral region than a single-object survey at a similar res¬ 
olution, on the order of 100 A for a multi-object (e.g., 
iFi'iresz e t al.||2008ll versus 1500 A for a single object sur- 
vev ('e.g.. lMarcv fc Bntleill199611 . This has the immediate 
implication that choosing the wavelength range used by 
a multi-object survey will be much more important than 
in a single-object survey. Secondly, a multi-object RV 
survey may not conduct an extensive pre-survey vetting 
of targets. This will give their target lists a higher dis¬ 
persion in stellar mass, rotation velocities, activity, sur¬ 
face gravity, and metallicity than current single-object 
surveys. Finally, multi-object surveys will operate in an 
observing mode more similar to photometric transit sur¬ 
veys for exoplanets, as compared to traditional RV sur¬ 
veys. Instead of exposing to reach the same SNR for 
each of the survey’s target stars, a multi-object survey 
will expose for the same time on each star. This will 
make the detection sensitivities of a multi-object funda¬ 
mentally different from a single-object survey in terms 
of stellar mass, effective temperature, and metallicity. A 
proper understanding of how these sensitivities change is 
vital to the initial design of a multi-object survey, and 
understanding the exoplanet statistics from the survey 
when it has completed. 

There have been several efforts to study the sources of 
radial velocity uncertainty since the beginnings of pre¬ 
cision RV surveys. The first description of how to nu¬ 
merically calculate the expected velocity uncertainty us- 
ing the p roper ties of an observed spectrum was given by 
iConnesI (II985I1 . and who also introduced the “Q” fac¬ 
tor as a way of parameterizing the amount of RV in¬ 
formation available in a spectral region, though they do 
not discuss how this is expected to vary with spectral 
properties. [Murdoch fc Hearnshawl (|I99If) performed a 
similar analysis, but started from the idea of measuring 
cross-correlatio n functions, ra ther than measuring spec¬ 
tral slopes as in lConn^ (I1985D . which makes for an inter¬ 
esting conceptual comp arison b e tween these two works. 
The approach used by IConnesI was repeated in 

iButler et al.l (|1996ll . who were interested in comparing 
theoretical Poisson-limited velocity unce rtainties agains t 
their actual observations. Similar to IConnesI (II985I1 . 
they did not discu ss how the u n certai nty varies with 
stellar properties. iBouchv et al.l (120011) . after present¬ 
ing their own derivation of the photon-limited RV preci¬ 
sion, were the first to examine how velocity uncertainty 
changes as a function of spectral type, stellar rotation 
velocity, and spectroscopic resolution. Again, these re¬ 
sults were used as part of a g eneral discussion o f the 
capabilities of the CORAL IE (iQuel oz et al. l I2000D and 
the (the n future) HARPS (IWavor et al.l 120031) ^ectro- 
graphs. IBouchv et al.l (1200 ID thus concentrated their 
analysis on a limited range of spectral types, rotations. 


and resolutions, and drew descriptive conclusions from 
their results (for example, they inferred that uncertainty 
is prop ortional to v si n i whe n v sin i is large). More re¬ 
cently, [Bottom et al.l (|2013D considered the sources of 
uncertainty in RV measurements as part of an exami¬ 
nation into how to optimiz e RV s urveys of GKM stars . 
Unlike iButler et ah (|1996D and IBouchv et al.l (1200ID . 
[Bottom et al.l ( 2013 ) did not derive an equation for the 
photon-limited uncertainty in an RV observation, but 
took model spectra, added Gaussian noise, and then fit 
for the Doppler shift in the noisy spectra using cross¬ 
correlation techniques. The authors intent was to bet¬ 
ter replicate the true process b e hind t he measurement 
of stellar RVs. [Bottom et al.l (|2013D consider wider 


and more finely-spaced ranges of wavelen gth, tempera¬ 
ture, and spectral resolution compared to IBouchv et al.l 
(j200lD . but they do not consider in de tail the effect 
of cha n ging stellar ro t ation. Similarly to IBouchv et al.l 
(|200lD . [Bottom et al.l (I2013D restrict themselves to de¬ 
scriptive conclusions regarding the sources of velocity un¬ 
certainty. 

In this paper, we aim to provide a more thorough de¬ 
scription of photon-limited stellar velocity uncertainties. 
We do not consider the effects of “jitter” sources such as 
spots, granulation, or asteroseismic pulsations. Similarly, 
we do not treat instrumental velocity uncertainties like 
wavelength calibration, optical effects, or instrumental 
drifts, nor do we consider the effect of telluric absorption 
lines on the spectra. Instead, we are here focused on the 
statistical velocity uncertainty one would measure for a 
photospherically stable rotating star using a perfect in¬ 
strument through a completely transparent atmosphere. 

Starting from a basic derivation of how to calculate 
the velocity uncertainty in a spectrum, our intent is to 
consider, in detail, the effects of effective temperature, 
surface gravity, metallicity, stellar rotation, spectral reso¬ 
lution, and macroturbulence on the uncertainty in an RV 
measurement. This allows us to make not only descrip¬ 
tive, but prescriptive conclusions regarding the sources 
of uncertainty. For example, we discuss precisely why 
stellar rotation has the effect it does on velocity uncer¬ 
tainties, starting from the shape of the rotation kernel 
itself. We are thus able to simply and numerically de¬ 
scribe how rotation affects the velocity uncertainty across 
all rotation velocities. 

Our ultimate goal is to provide a simple and trans¬ 
parent description of the various sources of velocity un¬ 
certainty in RV measurements, so as to give the reader 
a clear picture of the interlocking forces at work. To 
this end, we use our results to provide several simplify¬ 
ing approximations that capture the dominant sources 
of velocity uncertainty as function of stellar mass. This 
allows us to better understand why current RV surveys 
achieve the precision they do, and to provide guidance 
for the design of future surveys. 


2. PHOTON-LIMITED RADIAL VELOCITY 
PRECISION 

Consider a one-dimensional spectrum A that is dis¬ 
cretely sampled at equal velocity intervals (e.g., bins or 
pixels) AV centered at velocities Vi. The number of pho¬ 
tons in each velocity bin AV is W = AAV. We will 
assume for simplicity that the spectrum is in units of 
photons per pixel or per resolution element, and the un- 
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certainties are Poisson-dominated so that the uncertainty 
in Ni is just (Ttv.i = We will further assume that 

these uncertainties are uncorrelated. 

Although the precise details vary depending on the 
method and instrumental setup, we can think of the pro¬ 
cess of determining the radial velocity of this spectrum 
as a cross-correlation of a template spectrum of the star 
(either from models or a previous observation) against 
the observed spectrum. Typically, the peak or centroid 
of the cross-correlation function output by the spectral 
template matching is assumed to be the best estimate of 
radial velocity of the spectrum relative to the template. 
We wish to determine what the uncertainty is in this es¬ 
timate. As described in the Introduction, this general 
problem has been considered before, and the following 
derivation of the general answer is similar to previous 
work. 

As mentioned, we will assume that the uncertainties 
in our imagined spectrum are purely from Poisson noise, 
are free of any systematics, and uncorrelated. Further¬ 
more, we will assume that the velocity offset is not co¬ 
variant with any other possible fitting parameters. We 
are therefore, in some sense, determining the best possi¬ 
ble precision on the velocity offset that can be obtained, 
given the structure of the spectrum and assuming only 
source photon noise. There are various ways to compute 
the lowest possible uncertainty i n the velocit y offset (i.e., 
the minimum variance bound, IGouldl [199511 . for exam¬ 
ple using a Fisher matrix estimation. However, we will 
simply use the maximum likelihood estimator and error 
propagation, which is mathematically equivalent to the 
Fisher estimate under our set of assumptions. 

We will assume that each discrete sampling of the spec¬ 
trum yields an estimate of the velocity Vi of the spectrum 
with an uncertainty and the best estimate of the ve¬ 
locity V of the entire spectrum is the weighted mean over 
all the samples, 




( 1 ) 


where Wi = 1/(7^. The variance in this mean is then 
determined from the general error propagation equation, 
i.e., by differentiating the maximum-likelihood function 
and keeping only the first order terms, such that. 


(T 


2 

V 


= a! 





( 2 ) 


From Equation (HD the partial derivatives in Equation 
m are then 


dV ^ Wj ^ 1/af 

Substituting and simplifying, we obtain (see also 
iBiitler et al.lll99(il l. 

4 = (V^) ^ 0) 

We can relate the uncertainty in the velocity inferred 
for each point in the spectrum ai to the uncertainty in 


the intensity (photon number) at that point cr/^i via the 
local derivative of the spectrum with velocity, 

_ a'N,t 

~ {dN/dV)U' 

Again assuming through Poisson statistics that = 
y/Nl, and subsisting the above into Equation (jlD 

. i 

This is generally how well we can measure the velocity 
positions of the features in an arbitrary function of in¬ 
tensity versus velocity. This is similar to the “Q” factor 
formulation first used stated in IConnesI (1198511 . 



2.1. Centroiding absorption lines 

Let us first consider the case of a Gaussian absorption 
line. For compactness, let us define G(V),Vbj0G) as the 
appropriately normalized Gaussian distribution centered 
at Vb and with a full-width at half-max (FWHM) 0^: 


G{v,Vo,eG) 



.0y(41n2)_ • 


(7) 


If we take a spectrum composed of points separated by 
a constant AU in velocity, and if this line absorbs Ntot 
photons, then a spectrum containing only this line can 
be described by 

N^iV) = (/o - NtotG[V,Vo, 0g])A1/, (8) 

where Jq is the continuum level, in units of photons per 
unit velocity, and NjiVi) is the number of photons in 
a particular velocity bin centered at V). The factor of 
AU is the velocity span of a pixel. We can rewrite this 
equation in terms of the velocity equivalent width of the 
line, W = Ntot/Io, as 

N^iV) = /oAU(l - WG) = A^,cont(l - WG). (9) 
Substituting into Equation ([6D, we get 


CTv = 




f 2 V 2 ln2\ 

^ I^AV^W^G^ 

1 j 

/oAU(l - lUG) 



. (10) 

By the Euler-Maclaurin formula, we can approximate 
this sum with the integral 


I 


/ OO 

(V-Vof 

-OO 


/ 2V2A^\ 


(11) 


I^AV^W^G^ ; 

IoAV{l-WG)^ 


Unfortunately, this integral has no analytic solution. As 
a limiting case, consider a shallow absorption line such 
that 1 — VUG ~ 1. Now we may analytically solve the 
above equation to get 


(7v = 


loW^ 

f 2 V 2 In 2^ 

4v^ 

^ 0G jj 


( 12 ) 
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Fig. 1. — The exponent with which the velocity uncertainty cry 
scales with the width © of an absorption line, i.e., cry oc 0^, 
as a function of the line depth relative to the continuum level. 
This is numerically calculated using Equation I© for each of the 
three line profiles. At a fractional depth of 5^ 0 we find that the 
uncertainty scales as ay oc 0^'^, which is the same as for ou r 
ana lytic a ppro ximations for shallow lines given in Equations lO, 
and GHii. 


or, 


0? 


\ f,3/2 

' '-’g 




0 


(v2(21n2)3/4 j ~ 0-69 vr/Fo' 


(13) 


In numerical tests (Figure [1}, we find that this approxi¬ 
mation is valid for lines with depths less than about 10% 
of the continuum level. For lines deeper than this the 
exponential dependence on the width of the Gaussian 
increases. By the time the depth of the line is 95% of 
the continuum level we find that the uncertainty in the 
centroid scales roughly as 0 q. 

The uncertainty in the centroid of a shallow Gaus¬ 
sian absorption line therefore scales approximately as 
the FWHM of the l ine to the three-halves (see also 
iLovis fc Fischeii 120111) . We can understand this intu¬ 
itively by approximating an absorption line as triangular 
in shape, with a peak that is lower than the continuum 
by a fractional amount 1/0, and a width of 0. The form 
for such a line is simply: 


N^{Vi)=IoAV 


(^1-W 


V, - Vo-e 

02 


(14) 


for \Vi — Vo| < 0 and N^{Vi) = IqAV otherwise. For 
a weak line, the second term in parenthesis in Equation 
(d is negligible, and applying Equation ([6]) yields 


1 03/^ 03/^ 

av = —7=- ~ 0.71-=, 

s/2Wy/To W/Fo 


(15) 


with the same scaling and very similar coefficient as the 
weak Gaussian absorption line case. Conceptually, we 


can understand this as follows. The amount of infor¬ 
mation decreases, and the uncertainty increases, as the 
inverse of slope of the triangular line: uy oc 0^. On the 
other hand, the amount of information increases, and 
the uncertainty decreases, as the square-root of the to¬ 
tal number of photons. In the approximation of a weak 
line this is proportional to the square-root of the velocity 
span of the continuum covered by the base of the trian¬ 
gular line, and thus cry oc 0“3/2^ xhe net result of these 
two competing effects is that cry oc 03/^. This general 
argument implies that uncertainty in the centroid for all 
weak lines in a background continuum, regardless of their 
detailed prohle, should scale as 03/^. We demonstrate 
this explicitly for the other velocity profiles we consider 
below. We note that this scaling contrasts with the un¬ 
certainty of measuring the centroid of a pure Gaussian 
profile with zero continuum (i.e., calculating the stan- 
dard error of the rnean), which has the scaling of (e.g., 
iKennev fc Keeping 1 19471) 


g _ 1 Q 

2y2h/2VW’ 

where a is the standard deviation of the Gaussian profile, 
and n is the number of observations in the sample. 

Similarly to the Gaussian, we can work through the 
corresponding derivation for a Lorentzian absorption line 
with FWHM 0i to find 


CTy 



1.25 


TF/Fo' 


(17) 


This analytic solution is, again, under the assumption 
that the line depth is negligible relative to the contin¬ 
uum level. As can be seen in Figure [1] the exponent on 
Ql increases as the line depth increases, similar to the 
Gaussian case. 

Often, spectral lines are effectively described with 
Voigt profiles: the convolution of a Lorentzian and a 
Gaussian. Though there is no analytic description of 
a Voigt profile, we numerically examined how the uncer- 
taint y in the centroid scales with t he Voigt width. We 
used lOlivero fc Longbothuml (|1977[) ’s approximation for 
the effective FWHM of a Voigt profile, 


0voigt = O.53460L + y^O.21660| + 0^, (18) 


where 0 l is the FWHM of the Lorentzian component 
and 0G is the Gaussian FWHM. We directly calculated 
the uncertainty in the centroid using Equation ([6]) for 
various values of 0 l and Qq. As with the pure Gaus¬ 
sian and pure Lorentzian, we find that the uncertainty in 
measuring the centroid of a Voigt profile is proportional 
to 0voigf For the constant of proportionality relating 
CTy and 0Voigt, we found numerically that it varied as a 
function of the ratio 0l/0g- 


cry = 0.96 



03/2 

'^Voigt 

wVh' 


(19) 


This form of the leading constant is good to 10% over 
the range of 1/5 < 0 l/0g < 5. 

In addition to spectral absorption lines, we will also 
need to consider the role stellar rotation plays in setting 
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measured RV uncertainties. We begin by considering the 
shape of the kernel itself to determine the FWHM of the 
rotation kernel: Orot- To do so, we must first assume 
a limb-darkening law. For simplicity, we use a simple 
linear limb-darkening law: / = /o(l — e-l- ecosff), where 
Iq is the intensity at the center of the stellar disk, 0 is 
the angle of the surface to our line of sight, and e is the 
limb-darkening coefficient. Following iGravl (|2008ll . the 
normalized rotation kernel is then 


G(Au) = 


2(1 - e)v^l - {Av/vr^ 
7TO™t(l - e/3) 
^ 7 re(l - {Av/vrotf ) 
7rurot(l - e/3) 


( 20 ) 


where Vrot is the rotation speed at the limb of the star. 
For the case of no limb-darkening (e = 0) this reduces to 


G(Au) =-\/l - (An/ur-ot)^- 


'^Vrot 


( 21 ) 


To find ©rot we then set G{Av) = l/iWrot (i-e., half the 
maximum), Av = 0rot/2, and solve. Thus 


©rot = Vs Vrot (for e = 0). (22) 

At the other extreme of e = 1, we may solve Equation 
(POl) to find 


©rot = V2 Vrot (for e = 1). (23) 

Aside for the cases of e = 0 and e = 1, Equation (1^ 
allows for no simple analytic formula for ©rot as a func¬ 
tion of Vrot and e. We therefore numerically measured the 
FWHM of several calculated kernels between 0 < e < 1. 
We found that in between the two limb-darkening ex¬ 
tremes the FWHM went linearly with e, such that 


©rot = [{V2-Vs)e + ^3] Urot. (24) 

This relation is accurate to better than 5% over the en¬ 
tire range of 0 < e < 1. For reference, a Sun-like star 
observed at 5500 A would have e « 0.75, and e ~ 0.4 if 
observed at 10,000 A. 

Now, similar to the absorption line profiles, we may 
use Equation ([6l) to determine how velocity uncertainty 
scales with the width of the rotation kernel. We first 
consider the case of a fully limb-darkened kernel with 
e = 1 and equivalent width W subtracted from a con¬ 
tinuum. Thus the spectrum is given by Nj{Vi) = 
A^7,cont(l — WG). This represents the ideal case of a 5- 
function absorption line being rotationally broadened by 
the kernel. After making the appropriate substitutions 
into Equation ([6]), and again assuming that \ — WG~ 1, 
we find that 


av = 


2/3 ©^/? 




23/2 


0.49 


© 


WVh 


(for e = 1). (25) 


Unfortunately, e = 1 is the only case for which we 
may calculate ay directly from the rotation kernel itself 
using Equation ®. For all other values of e the slope of 
the kernel goes to infinity as AU — Vrot- We were not 
able to find an appropriate analytic or numeric integral 
to avoid this, so we instead convolved kernels for e < 1 
with a normalized Gaussian of small fixed width (cr = 0.1 



Fig. 2.— The colored points show how the velocity uncer¬ 
tainty varies with FWHM for absorption lines with a Gaussian, 
Lorentzian, Voigt, or rotational profile with W = 0.2 km s“^ and 
7o = 1 photon per km s~^, determined via direct numerical com¬ 
putation using Equation The overplotted lines are the directly 
calculated velocity uncertainties for these same abs o rpti on lin es us¬ 
ing our analytic approximations in Equations lO, inii, lO, and 
llS}. For the Voigt profile we set the Lorenz tian and Gaussian 
widths equal, so that ©l/©g = 1 in Equation (|19ll . 


km s 3) and measured the velocity uncertainties of the 
resulting lines for Vrot > 25 km s”^. We found that 
the velocity uncertainty continued to be proportional to 

Orgt ) with the constant of proportionality varying with 
roughly linearly e, such that 


©3/2 

av^ (0.347+ 0.U6e)-^ 

WyGo 


(26) 


This is accurate to 2% over 0 < e < 1. 

Interestingly, if one rewrites Equation (j26p in terms of 
Vrot, rather than Orot, the e dependence of ay nearly can¬ 
cels out. Put another way, velocity uncertainties are not 
strongly effected by the precise amount of limb-darkening 
in the stellar photosphere. The difference in ay between 
e = 0 and e = 0.75 (the locations of the minimum and 
maximum of the proportionality coefficient) is only 5%. 
This is somewhat dependent on our choice of a linear 
limb-darkening law, but this result should approximately 
hold for more complicated limb-darkening laws. 

Figure [5] shows the directly calculated uncertainties for 
a rotation kernel, Gaussian, Lorentzian, and Voigt profile 
with the same equivalent width as a function of FWHM. 
The overplotted lines are what we expect for the un¬ 
certainty based on Equations (USD, (EZl), and (HU). For 
reference, Figure [3] shows all four profiles, each one with 
© = 1 km s“^ and W = 0.2 km s“^. 


3. GALGULATING STELLAR VELOGITY 
UNGERTAINTIES 
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Our ultimate goal is to calculate the photon-limited 
uncertainty expected for RV observations of a main se¬ 
quence star with a certain exposure time. To achieve 
this, we will examine how the velocity uncertainties of 
synthetic spectra - calculated directly using Equation 
m - change as a function of instrumental and stellar 
properties. So far we have been describing the behavior 
of individual lines; going forward we will instead consider 
the average behavior of groups of lines over an extended 
spectrum. 

From our consideration of the uncertainty in centroid- 
ing single absorption lines, we anticipate that the velocity 
uncertainty of a group of lines should scale as 


cry oc 


^/N^ 


lines 


03/2 


(27) 


where Iq is the continuum flux level, 0 and W are the 
average FWHM and equivalent widths of the lines in the 
spectrum of interest, and Mines is the number of lines in 
the spectrum. The continuum level is set by the spectral 
energy distribution of the target star, the exposure time, 
the size of the telescope, and the throughput of the opti¬ 
cal system. The average FWHM of the spectral lines is 
determined by the spectral resolution of the instrument, 
the observed wavelength range, as well as a host of stel¬ 
lar parameters such as mass, temperature, age, rotational 
velocity and metallicity. We chose to restrict ourselves to 
stars that are on the main sequence, specifically dwarfs 
from 2600K to 7600K. This allows us to parameterize 
stellar properties in terms of one parameter (e.g., mass 
or effective temperature). This leaves us with seven main 
parameters that will determine the average FWHM of 
the spectral lines: spectral resolution, wavelength range, 
stellar rotation, photospheric macroturbulence, effective 
temperature, surface gravity, and metallicity. 



Av (km s“') 

Fig. 3.— The different line broadening profiles we consider in 
Section 2.1, all with the same FWHM (0 = 1km s“^) and equiv¬ 
alent width (W = 0.2km s“^). 


For simplicity, we chose to quantify Mines, the num¬ 
ber of lines present in a spectrum, and W, their average 
equivalent widths, using separate and purely descriptive 
functions of effective temperature, surface gravity, and 
metallicity. Thus Equation (l?7l) becomes 

ay cx ^ filogg) /([Fe/H]). (28) 

v-'o 

This is the approximation we will use in our fitting. 

We are making an important assumption here: that 0 
exists as a meaningful descriptor for a given spectrum. 
We intend to calculate 0 for a spectrum as a combina¬ 
tion of some 00 caused by the inherent widths of the 
stellar lines, and an additional line width caused by one 
or more broadening mechanisms (e.g., Qr for spectral 
resolution). The assumption we are making therefore al¬ 
lows us to determine 0o not by measuring the FWHM 
of each line in a given spectrum, but rather by fitting 
how that spectrum’s velocity precision scales as a source 
of broadening is applied to it. If this assumption that 
0 exists as a meaningful quantity is correct, then we ex¬ 
pect that ay for spectra of varying wavelengths will be¬ 
have roughly self-similarly as we apply various broaden¬ 
ing mechanisms, and that spectra over a wide wavelength 
range will have fitted values of 0o that are similar. For 
the moment, let us accept this underlying assumption as 
correct. We shall see in subsequent subsections that ay 
does indeed behave self-similarly, and that a single value 
of 00 can describe a spectrum’s response to broadening 
for different spectral regions spaced over thousands of 
angstroms. 

We used two different sets of synthetic spectra in 
our fitting: one s et from the BT-Settl model spectra 
(|Allard et al.ll2012ll and another ca lculated using Kurucz 
model atmospheres (|Kurucdll992L hereafter Kurucz92). 
The two sets provide us with different pieces of infor¬ 
mation: the BT-Settl spectra cover the full temperature 
range of interest, but have a (relatively) coarse wave¬ 
length spacing, while the Kurucz92 spectra are extremely 
finely spaced in wavelength but are only available for a 
subset of the temperatures in which are interested. We 
therefore used the BT-Settl spectra to examine the ef¬ 
fect of stellar effective temperature and surface gravity, 
and used the Kurucz92 spectra to model line broaden¬ 
ing mechanisms like spectral resolution and stellar ro¬ 
tation. Since one can generically think of these latter 
mechanisms as externalities imposed upon “perfect” stel¬ 
lar spectra (i.e., rotational broadening is not an intrinsic 
part of creating absorption lines) we expected the results 
we find using the Kurucz92 spectra to be consistent with 
our results using the BT-Settl spectra, once we correct 
for the difference in wavelength sampling. As described 
later, we ultimately found this to be the case. 

We used flux-normalized spectra in three broad bands: 
“optical” spectra from 4000 A to 6500 A, “red” spec¬ 
tra from 6500 A to 10000 A, and “near-infrared” (NIR) 
spectra from 10000 A to 25000 A. The BT-Settl spectra 
used the lAsplund et ahl (|2009f ) solar abundances, had so¬ 
lar metallicity with no a-enhancement, and were spaced 
200K apart from 2600K to 7600K. The BT-Settle spectra 
available for download use a variable wavelength spac¬ 
ing, with a finer spacing occurring around the absorp¬ 
tion lines. On average the wavelength spacing was 0.05 A 
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(i? Ri 100,000) in the optical, 0.05 A {R fa 165,000) in 
the red, and 0.2 A (i? fa 88, 000) in the NIR. 

For the Kurucz92 spectra we used the ODFNEW ver¬ 
sions of the Kurucz92 model atmospheres with no a- 
enhancement, and gene rated the spectra with v2.76 of 
iGrav fc Corballvl (|1994f) ’s spectrum code. For the op¬ 
tical and red spectra we used a fixed wavelength spacing 
of 0.001 A {R fa 6 X 10®) and for the NIR spectra we 
used 0.005 A (i? Ri 3.5 x 10®). In all the bands we set the 
microturbulent velocity to 1 km s“^ and left the macro- 
turbulent velocity at zero. We considered the effect of 
macroturbulence separately. The Kurucz92 models cov¬ 
ered effective temperatures from 4000K to 7500K with a 
spacing of 250K. 

To investigate wavelength dependent features, we di¬ 
vided all our spectra into 100 A chunks. This partially 
isolates individual line groupings, like the Mg B triplet, 
so that we can test whether these groupings react to 
changes in a self-similar way. Furthermore, splitting the 
spectra into 100 A chunks allowed us to mirror the actual 
analysis procedures of current multi-order RV surveys, 
and is representative of the amount of spectral informa¬ 
tion available in proposed multi-object surveys. 

For each chunk, we calculated the expected velocity 
uncertainty using Equation (jT]) . After transforming each 
chunk from wavelength to velocity space, we normalized 
the chunks so that each had Ni = 1 photon per m s“^ in 
the continuum (/q = 1). We then numerically calculated 
the slope of the spectrum at each pixel {{dN/dV)\i). By 
Equation (|3]), this then gives us a velocity uncertainty 
for each 100 A chunk. 

We now wish to see how the velocity uncertainties cal¬ 
culated from the individual 100 A chunks vary as we vary 
spectral resolution, stellar rotation, macroturbulence, ef¬ 
fective temperature, surface gravity, and metallicity. 

3.1. Spectral Resolution, Stellar Rotation, and 
Macroturbulence: Kurucz92 Based 

Conceptually, the effect on a spectrum of changing 
spectral resolution (R), stellar rotation (?;sini) or the 
macroturbulent velocity {vmac) can be viewed as an ex¬ 
ternality imposed upon a “perfect” spectrum with R = 
oo, wsinj=0 and Vmac=0- Regardless of the underlying 
stellar parameters, the velocity uncertainty of a spectrum 
should vary (roughly) with the same functional form for 
R, tisini and Vmac- We therefore rewrite Equation (I27|) 
to 

'^mac)] ^ 

where ctv.o is the velocity uncertainty of a “perfect” spec¬ 
trum with R=oo, and no rotation or macroturbulence. 
(preiiR,Vrot,Vmac) IS the increase in the average FWHM 
of the spectral lines caused by changes in R, z;sin*, and 
macroturbulence relative to that of the perfect spec¬ 
trum. We defined iprei such that 0; 0)=1 and 

(/JreAO; oo, 0) = (/JreAO; 0, oo) = OO. We assumed that the 
R, Vrot and Vmac contributions to iprei were separable, 
and we find that this is approximately true. 

We first considered the R dependence of (prei using a 
Kurucz92-based spectrum of a a 5750K, log( 5 ) = 4.5, 
[Ee/H]=0.0, Sun-like star split up into 100 A spectral 
chunks. We used a Kurucz92-based spectrum - instead 
of a BT-Settl spectrum - because of the extremely fine 


! Optical (4000-65001) ! 

Red (6500-10000l) 

; NIE (10000-250001) 1 








^ 10 '* 10 -“^ 10 " 10 -® 10 " 10 ® 

Spectral Resolution (R) 


Fig. 4.— The velocity uncertainty of 100 A chunks as a function 
of spectral resolution using Kurucz92-based spectra. The top row 
shows the absolute uncertainty for the chunks in our three wave¬ 
length regions, while the middle row shows the velocity uncertainty 
for each chunk normalized to /?=50,000. The red line in the middle 
panel is the median of the relative velocity uncertainties. The bot¬ 
tom row shows the fractional difference between each individual 
chunk’s relative uncertainties and this median. The middle row 
demonstrates that all of the wavelength chunks respond roughly 
self-similarly to changes in spectral resolution. 


wavelength spacing available with the Kurucz92 spec¬ 
tra. Our Kurucz92 spectrum had a wavelength spacing 
of 0.001 A, as compared to a median spacing of 0.05 A 
in the BT-Settl spectrum. Though the BT-Settl spacing 
gives a well sampled spectrum for most applications, we 
will see that for our specific examination of line broad¬ 
ening mechanisms the 0.05 A spacing has a noticeable 
effect, by effectively setting a base spectral resolution of 
i?Ril05,000 (see Figure[5|). 

We assumed that the effect of instrumental spectral 
resolution could be approximated by convolving a spec¬ 
trum with a Gaussian of EWHM equal to c/R. We there¬ 
fore convolved each 100 A chunk with Gaussians corre¬ 
sponding to a range of i?-values, and calculated the ve¬ 
locity uncertainties. 

To do this, we first chose a continuum SNR per ve¬ 
locity bin of unity, by setting = 1 in the continuum 
(and with corresponding lower values of Ni in the ab¬ 
sorbed portions of the spectrum). We then applied a 
given spectral resolution {(prei[R, 0,0]) via Gaussian con¬ 
volution, calculated the slope at each bin in the synthetic 
spectrum, {dN/dV)\i, and used Equation (g]) to calculate 
the velocity uncertainty for the chunk as a function of R. 

Eigure HI shows the absolute values of the chunk ve¬ 
locity uncertainties for a 5750K star from i?=10,000 to 
i?=200,000 in the top row. We have separated the chunks 
into “optical”, “red” and “NIR”. In absolute terms there 
is a range in uncertainties across the chunks as a result 
of specific spectral features in specific locations. 

To check for any wavelength dependencies, we normal- 
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ized the curves from the top row of Figure|3]by their value 
at _R=50,000. The middle row of Figure |4] shows these 
normalized curves, and the red line shows their median 
value at each spectral resolution. One effect to note is 
that the dispersion of the chunks in the middle row pan¬ 
els increases as one moves to the red; that is, while the 
optical chunks all behave very similarly, the NIR chunks 
show more variation relative to each other. Specifically, 
many of the chunks in the NIR seem to be less affected 
by changing resolution than the optical chunks. 

This occurs because the less affected chunks, which are 
mostly at longer wavelengths, have residual molecular 
features in the model spectra we are using. Since the lines 
that make up the molecular bandheads are very closely 
spaced, they blend together into composite lines with a 
large widths even at high resolutions. This lessens the 
effect of increasing spectral resolution in resolving these 
features. 

The red line in the middle row of Figure|3]is the median 
relative velocity uncertainty across all the chunks, and 
the bottom row of Figure|4]shows the fractional difference 
between all the chunks and this median. The fractional 
difference across the entire optical and red wavelength 
ranges is rarely more than 20%, while the NIR chunks 
stay within about 30% of the calculated chunk median. 

We repeated the above procedure to numerically calcu¬ 
late the median velocity uncertainty vs. R using all the 
Kurucz92-derived spectra from 4000K to 7500K in steps 
of 250K. Across this temperature range the results for 
the optical, red, and NIR chunks were similar to our il¬ 
lustrative, 5750K, example. The velocity uncertainties of 
the chunks roughly stayed within 25% of the calculated 
median. 

To determine the R dependence of (prei, we fit to the 
chunk median velocity uncertainties. Since we are inter¬ 
ested in the relative change in velocity uncertainty for the 
chunks, we re-normalized each chunk median so that the 
velocity uncertainty at i?=3 x 10® was unity, such that 
(frei{oo, 0,0) = 1. To describe the average FWHM of the 
spectral lines in the chunks, we fit the chunk median as 
a Voigt profile with some inherent width 0o, such that 
the relative increase in the average FWHM of the chunks 

scaled as with 

O.53460O + v'0.21660^ + 0^ . 

T’r-ei — „ I (3Uj 

<=>0 

where Qr = c/R. For each wavelength region, we fit the 
measured median velocity uncertainties as a function of 
resolution using Equations (1^ and (|5(I)) and by finding 
the best-fit value of 00- Figure[5]shows the median chunk 
velocity uncertainty in the optical for a 5750K star as a 
function of Qr in black, overlaid by with our best-fit in 
green. The residuals to the best-fit are no more than 3% 
across the entire range of resolutions at 5750K. 

As an illustration of the effects of the different wave¬ 
length spacing in the Kurucz92 spectra and the BT-Settl 
spectra. Figure [5] also shows a similarly calculated curve 
for the velocity uncertainty of a 5800K, log(g) = 4.5 
BT-Settl spectrum in red. Note that the BT-Settl curve 
asymptotes to a significantly higher velocity uncertainty 
at high resolution (low 0 _r). We interpret this as a 
result of the coarser wavelength spacing in the BT- 
Settl spectrum, which imposes a base “resolution” of 



Fig. 5.— Velocity uncertainty, relative to the velocity uncertainty 
in a Kurucz92-based spectra at R=S x 10®, plotted as a function 
of Qfi = c/R. The black points show how the uncertainty on 
Kurucz92-based spectra with a wavelength spacing of AA = 0.001 
varies with re solu tion, and the green line shows our best fit us¬ 
ing Equation 1301 . In addition, note the difference between the 
black Kurucz92-based spectra and a similar analysis done on the 
BT-Settl-based spectra (red points). The wavelength spacing of 
the BT-Settl spectra (AA ~ 0.05) imposes a resolution “floor” of 
/l?iil05,000. This is demonstrated by the orange points, which 
show the behavior of the BT-Settl spectra if if we include a base of 
il=105,000 and add this in quadrature to the resolution increase 
being applied by our Gaussian convolution, thereby imposing a res¬ 
olution floor of ©/? = 2.9 km s“^. This transformed BT-Settl be¬ 
haves similarly to the black Kurucz92 line. Alternatively, the blue 
points show how a Kurucz92 spectrum sampled with a AA = 0.05 
wavelength spacing, comparable to the BT-Settl spectra, reacts to 
changing spectral resolution. 


R Ri 5250A/O.O5A = 105,000 in the optical. Indeed, if 
we include a base of 77=105,000 and add this in quadra¬ 
ture to the resolution increase being applied by our Gaus¬ 
sian convolution, thereby imposing a resolution floor of 
Qr = 2.9km s“^, the red BT-Settl curve in Figure [5] 
transforms to the orange line in Figure [5] and nearly 
matches the Kurucz92-based results. As a further test 
we also generated a Kurucz92 spectrum with a 0.05 A 
wavelength spacing, and its curve is shown in blue. This 
coarser Kurucz92 spectrum nearly matches the BT-Settl 
results, which makes us confident that the difference be¬ 
tween the 0.001 A Kurucz92 spectrum and the 0.05 A 
BT-Settl spectrum is primarily a result of the different 
wavelength spacings. Figure [5] also illustrates why we 
used the Kurucz92-based spectra for our examination 
of line broadening mechanisms: the unbroadened BT- 
Settl spectra are not sampled finely enough to represent 

V5reAoO,0, 0) = 1. 

We note that at low-77 in Figure [5] there is an offset 
between the Kurucz92 and BT-Settl spectra. We were 
not able to completely determine the cause of this offset, 
which is equal to about 0.15 dex at 77=3,000 and about 
0.05 dex at 77=30,000. 
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Fig. 6.— The median velocity uncertainty of a 5750K, log^ = 4.5 
star, relative to the velocity uncertainty for that same star at 
_R=lxlO®, as a function of the width of the underlying Voigt profile 
(Equation I18I V This figure highlights that the velocity uncertainty 
is not set by the resolution alone, but by the combination of the un¬ 
derlying average line width and the spectral resolution. The Voigt 
width used here is the combination of t he average inherent line 
width in the spectrum from Equation lEU and a changing spectral 
resolution, from R=1 x 10® to i?=10,000. While the calculated 
3/2 

points deviate from cry oc ^y^igt ^ small scale, the overall best 

fit to the median relative velocities is given by a ^yo%t ^ ^Voigt 
scaling. 


3.1.1. Comparison to Other Work 

It is also worth noting at this point that other au- 
thors’ (e.g..lHatzes fc CochraDlll99a iBouchv et al.ll2001l : 

iBottom et alJ I2013D numerical calculations of the de¬ 
pendence of how velocity uncertainty scales with spec¬ 
tral resolution find that at lo w-R the uncertainty goes 
approximately as ay oc R~^ (iHatzes fc Cochraiil 119921 : 
iBouchv et al.ll200Tll or ay oc R~^ ‘^ (iBottom et al.l[20^ . 
There are two connected points to consider here. First, 
as we have seen, we mathematically expect the velocity 
uncertainty to scale as 0^'^, where 0 is set by both Qji 
and the inherent line width 0o. This means that con¬ 
sidering only 0/i, as these authors do, does not account 
for the effect of the intrinsic width of the lines on the 
velocity uncertainty. 

Second, as we shall see, for a Sun-like star in the 
optical, 0/j dominates the inherent line width (i.e., 
0B lO0n) only fo r re solutions less then 6 ,000. Thus 
iBoudiv et al.l (1200111 and IBottom et al.l (|2013ll . who con¬ 
sider down to i?=10,000, find ay oc R~^ and ay oc R~^ ‘^, 
respectively, since they are largely fit ting over the tran¬ 
sition regime between On and ©o- IHatzes fc Cochraiil 
(jl992ll directly measure the velocity uncertainty at 
i?=2,500 and find ay oc R~^, but we believe this result to 
be a poor fit to their measurements, since this line passes 
substantially underneath the R=2,50 0 point. Indeed, if 
we take the three points in Figure 1 of lHatzes fc CochranI 



Optical 


4000 4500 5000 5500 6000 6500 7000 7500 


Teff (K) 


Fig. 7.— The ratio of our fits f or th e velocity uncertainties, calcu¬ 
lated using Equations ll30l l and l|3TJ, to the numerically computed 
chunk uncertainties as a function of changing spectral resolution 
and effective temperature in the optical. The RMS across all reso¬ 
lutions and temperatures is 7%, and the peak value is about 20%. 


(|1992ll and fit them using our formalism, we recover 
ay oc R~^'^ at low resolution. 

If we plot the median relative velocity uncertainty for 
a 5750K, log 5 = 4.5, star against the underlying Voigt 
width (see Equation [18] ) instead of the width solely due 
to spectral resolution (Figure (Bj), we can immediately see 
that the best fit to the calculated uncertainties goes as 
ay oc 0^/^. 

3.1.2. Temperature Dependence of Qo 

We repeated our procedure of fitting the velocity un¬ 
certainty as a function of Qr for the full tempera¬ 
ture range covered by the Kurucz92 models (4000K to 
7500K), and for the three wavelength ranges we consider 
(optical, red, and NIR). We found that the best-fit value 
of 00 decreased roughly linearly with temperature and 
was slightly different in each regime. Specifically, 0o goes 
as 

00 = 5.10521 km s"^ (1 - 0.6395 AT^ff) (31) 

for 4000 to 6500 A 

00 = 3.73956 km s"^ (1 - 0.1449 ATe//) 

for 6500 to 10000 A 

00 = 6.42622 km s"^ (1 - 0.2737 ATe//) 

for 10000 to 25000 A, 

where ATg// = Te///5800K— 1. These equations are 
good to 3% in the optical and 5% in the red and NIR. 

Figures mu and [HI show the ratio of our fits to the 
numerically calculated chunk velocity uncertainties as a 
function of effective temperature and spectral resolution. 
Across the entire range of temperature and resolution the 
ratios have a standard deviation of about 5% about unity. 
The largest differences occur in the optical at very low 



































10 


Beatty & Gaudi 



NIR 


4000 4500 5000 5500 6000 6500 7000 7500 


4.0 


Red 


4000 4500 5000 5500 6000 6500 7000 7500 

Teff (K) Teff (K) 



1.05 


Fig. 8.— The ratio of our fits f or th e velocity uncertainties, calcu¬ 
lated using Equations iIMt and eu, to the numerically computed 
chunk uncertainties as a function of changing spectral resolution 
and effective temperature in the optical. The RMS across all reso¬ 
lutions and temperatures is 4%, and the peak value is about 10%. 

and very high temperatures for resolutions near 60,000. 
Our fits under-predict the uncertainties at low temper¬ 
atures and over-predict at high temperatures, both by 
about 20% at these specific locations. 


3.1.3. Stellar Rotation 

We next turned to the effect of stellar rotation on the 
velocity uncertainty. We presumed that the effect of 
changing Vrot is similar to spectral resolution R, in that 
convolving a spectrum with a rotation kernel is similar 
to convolving a spectrum with a Gaussian, and that we 
could describe the change in the line width similarly as 




O.53460O + ^0.21660^ -b 

0 ^ 


(32) 


Here &rot is the FWHM of the rotation kernel. Equation 
(EH) presumes a Gaussian FWHM measurement, so for 
the purposes of Equation (1321) we will need to calculate 
an “equivalent Gaussian width” for the rotation kernel. 
This is the width of a Gaussian which, when convolved 
with a spectrum, will cause the equivalent broadening 
and associated increase in the velocity uncertainty as 
the rotational kernel. As we have seen in Section 2.1, 
this substitution may be easily made: the increase in ve¬ 
locity uncertainty caused by Gaussian broadening differs 
from rotational broadening only by a leading numerical 
coefficient. 

To determine the equivalent Gaussian width of the ro¬ 
tation kernel, we set Equations m and (12611 equal to 
each other and solve for 0^. Thus, 


0G,e 


0.347-b 0.146 
069 


0 ^ 


As we noted in Section 2.1, the velocity uncertainty 
caused by the rotational broadening at a particular ro- 


Fig. 9.— The ratio of our fits f or th e velocity uncertainties, calcu¬ 
lated using Equations iIMt and eu, to the numerically computed 
chunk uncertainties as a function of changing spectral resolution 
and effective temperature in the optical. The RMS across all res¬ 
olutions and temperatures is 3.5%, and the peak value is about 
10 %. 

tation velocity is only weakly dependent on the precise 
value of e. The difference in the proportionality constant 
relating 0%^ and ay varies by about 10% from minimum 
to maximum. We will therefore take the average value, 
which occurs at e = 0.5, for all of our results. This makes 
0G.eg = 0.72 Qrot- 

Having determined 0G.eij for rotation, we know wish 
to know how velocity uncertainty scales with stellar rota¬ 
tion velocity. In a manner similar to how we approached 
spectral resolution, we calculated how the velocity un¬ 
certainty in the wavelength chunks changed as Vrot went 
from 0 to 25 km/s using Kurucz92-based spectra. We 
used the AVSINI routine packaged with the spectrum 
code to apply the rotation kernel to our spectra using 
with e = 0.5. As one can see in Figure [101 changing Vrot 
is similar to changing spectral resolution in that it largely 
effects all of the wavelength chunks in the same way. The 
middle row of Figure [TUI shows the relative change in the 
velocity uncertainty normalized to Vrot = 5 km s“^ for a 
5750K, log( 5 ) = 4.5, [Fe/H]=0.0 star in our three wave¬ 
length regimes. Similar to our approach to fitting the 
effect of changing spectral resolution, we also calculated 
median values for our entire temperature range. 

We then fit the chunk averages in the same manner as 
for spectral resolution. In doing so, we found that rota¬ 
tional velocity affects the relative velocity uncertainty of 
a spectrum in almost exactly the same way as does spec¬ 
tral resolution. That is, when QG,eq = 0r the velocity 
uncertainty is almost exactly the same across our entire 
temperature range. 

3.1.4. Macroturbulence 

In addition to rotation, we also considered the effect 
of macroturbulence in the stellar atmosphere. For sim¬ 
plicity we assumed simple isotropic macroturbulence, so 
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that the effect of macroturbulence with velocity Vmac is 
the same as convolving a spectrum with a normalized 
Gaussian with standard deviation Vmacl‘^- Under this 
assumption the FWHM of the macroturbulence kernel is 
then simply 

e^nac = 2V^ ^ ^1.18 Vraac (34) 

Note that in reality, the effects of rotation and macro¬ 
turbulence are di fficult to observationally separate when 
Vrot ~ Vmac le.g.. IValenti fc Fischer! [200^ . This is par¬ 
tially a result of the fact that the effect of real macro¬ 
turbulence is not isotropic, and partially because macro- 
turbulent and rotational broadening are observed as a 
disk-integrated broadening profile. This makes the two 
effects difficult to separate observationally. 

Typical macroturbulent veloc ities for field dwarfs are 
on the ord e r of a few km s~^ (| Valenti fc Fisch^ 120051 : 
IGravl [2008t iBruntt et al.l I2010ll . with mid F-dwarfs at 
about 6 km s“^ and decreasing linearly with spectral 
type to about 1.5 km for an early K-dwarf. We used 
the empirical relatio n for Vm.ac as a func tion of temper¬ 
ature determined bv IBruntt et al.l (120101 1. which we list 
along with other stellar properties in Section 4. 

We therefore will use the same results we had for spec¬ 
tral resolution (Voigt line profiles, temperature depen¬ 
dence) and apply it to rotation and macroturbulence. 
Putting this all together, we may rewrite Equation (1291) 



1 10 1 10 1 10 


(km s"^) 

Fig. 10.— The velocity uncertainty of 100 A chunks as a func¬ 
tion of stellar rotation velocity from Kurucz92-based spectra. The 
top row shows the absolute uncertainty for the chunks in our three 
wavelength regions, while the middle row shows the velocity uncer¬ 
tainty for each chunk normalized to Vrot=^ km s~^. The red line 
in the middle panel is the median of the relative velocity uncer¬ 
tainties. The bottom row shows the fractional difference between 
each individual chunk’s relative uncertainties and this median. The 
middle row demonstrates that all of the wavelength chunks respond 
roughly self-similarly to changes in stellar rotation. 


as 


f O.53460o(re//) 

Iv 0o(Te//) 


(35) 


-(- 


v/0.21660^ + Q% + 0.518 8^,, + ^ 

QoiTeff) 


3/2 


3.2. Temperature: BT-Settl Based 

Temperature affects both the width and the number of 
lines usable for radial velocity measurements in a spec¬ 
trum, and unlike the line broadening mechanisms consid¬ 
ered in above, stellar temperature should be considered 
an intrinsic part of line generation. Without a detailed 
treatment of how spectral lines are created, it is there¬ 
fore difficult to arrive at a physically motivated analytic 
expression for how the velocity uncertainty in a spec¬ 
trum changes along with effective temperature. While 
the thermal velocity width of the lines will scale sim¬ 
ply as the square-root of the effective temperature, the 
pressure of the atmospheric layer where these lines are 
generated will change as well. These two competing ef¬ 
fects - temperature width and pressure width - are not 
easily separable. In addition, the number of lines in a 
spectrum depends upon a host of factors such as opaci¬ 
ties, atomic energy levels, and ionization equilibria that 
also provide no simple scaling with temperature. 

Figure Hi] is an illustration of the complexity of the is¬ 
sue. We selected three isolated iron lines, one silicon line, 
and one nickel line, and measured the FWHM of these 



Fig. 11.— Measured FWHMs of specific isolated lines as a func¬ 
tion of effective temperature. We have normalized each to be unity 
at 5800K, so as to identify any broad trends. The lines are color- 
coded according to the atom responsible: red are iron lines at 
5294.5 A, 5905.7 A, and 6078.5 A, blue is a silicon line at 6125.0 A, 
and green is a nickel line at 6482.7 A. While the FWHMs of all the 
lines generally grow smaller as temperature increases, one can see 
the large variation in how specific lines react. This is in contrast 
to, for example, spectral resolution, which will alter the FWHM of 
lines uniformly. 
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Fig. 12. — The velocity uncertainty of 100 A chunks as a function 
of stellar effective temperature using BT-Settl-based spectra. The 
top row shows the absolute uncertainty for the chunks in our three 
wavelength regions, while the middle row shows the velocity uncer¬ 
tainty for each chunk normalized that for 5800K. The red points 
in the middle panel is the median of the relative velocity uncer¬ 
tainties, while the orange overplotted line is our fit to the median. 
The bottom row shows the fractional difference between each indi¬ 
vidual chunk’s relative uncertainties and this median. The chunks 
show considerably more relative differences as a function of temper¬ 
ature compared to changing spectral resolution or rotation. This 
reflects the complicated competing effects that changing tempera¬ 
ture causes in spectral line generation, including changing pressure 
levels, ionization states, and opacities. 


lines as a function of temperature in a series of [Fe/H]=0, 
\og{g) = 4.5, Kurucz92 spectra. Our naive expectation 
was that the FWHM of all the lines would behave self- 
similarly and increase as the square-root of the effective 
temperature. Instead, the measured FWHMs decreased 
with effective temperature and display an idiosyncratic 
temperature dependence. As mentioned above, we at¬ 
tribute this divergence from our expected behavior to a 
set of competing effects, including varying local pressure, 
differences in ionization levels, and differences in opaci¬ 
ties. 

We therefore determined a purely numerical and de¬ 
scriptive scaling for how the relative velocity uncertainty 
in a spectrum changes with effective temperature. Due 
to their availability over a greater range of temperatures 
we used the BT-Settl models for this fitting. We again 
used 100 A chunks sliced out of spectra with effective 
temperatures of 2600K to 7600K and logg = 4.5 in the 
three wavelength ranges we are considering. 

The top panels of Figure IT^ show the effect of chang¬ 
ing temperature on the velocity uncertainty for each of 
the chunks. It is immediately apparent from Figure [12] 
that changing the effective temperature acts in a much 
less self-similar way across the chunks as compared to 
the external line-broadening mechanisms we considered 
previously. Not surprisingly, while the behavior of the 
optical chunks is roughly self-similar (left side of the mid¬ 


dle row of Figure [T^ . the NIR chunks (right side of the 
middle row of Figure [T^ show considerable differences. 
This is largely caused by the different line generation 
mechanisms at optical and NIR wavelengths. While the 
optical is mostly populated by atomic lines that change 
strength relatively slowly with effective temperature, the 
NIR chunks possess more molecular lines that have a 
sharp temperature dependence. For example, while a 
4000K, log((j') = 4.5, spectra from 24500 A to 24600 A is 
a forest of CO molecular lines, that same 100 A chunk in 
a 7000K, log(g) = 4.5, spectrum has only one atomic Fe 
and one atomic Mg line as its major spectral features. 

This also illustrates the vital importance of choos¬ 
ing the appropriate wavelength range in the NIR when 
designing an RV survey. For example, our illustrative 
24500 A to 24600 A chunk is a perfect example of a wave¬ 
length region that would be a reasonable choice for an 
RV survey focusing on K and M stars, but it would be 
a poor choice for a NIR survey that would observe FGK 
dwarfs. We consider the choice of wavelength range in 
more detail in the discussion section. 

To generally describe the behavior of the wavelength 
chunks as a function of temperature, we roughly approxi¬ 
mated the chunk medians by taking a least squares poly¬ 
nomial fit to the three wavelength regions, which yielded 

f{Teff)opt = 1 + 2.04515 ATe// + 3.13362 (36) 

+A.2S845AT^ff 

fiTeff)Red = 1 + 2.18311 ATe// + 4.00361 AT^ff 
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Fig. 13. — Velocity uncertainty as a function of effective temper¬ 
ature and surface gravity, relative to log( 5 ') = 4.5 using BT-Settl 
spectra. The black lines are individual 100 A spectral chunks, and 
the red points are the medians of all the chunks at the effective 
temperatures of the model spectra. The overplotted orange lines 
are our fits to the chunk medians. Notice that at the extremes 
of the temperate range we consider, where we expect the stellar 
surface gravity to be most different from log(p) = 4.5, the effect 
of changing surface gravity is the least. Surface gravity therefore 
plays a minor ('^10%) in setting velocity uncertainties. 
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+5.62077 

/(Te//)NIR = 1 + 1.62418 ATe// + 2.62018 AT^^ 

+5.01776 AT^ff, 

where ATg f f = Te///5800K — 1. The bottom panels 
of Figure [12] show the fractional difference between each 
chunk’s relative velocity uncertainties and the chunk me¬ 
dians. The relations in Equation (I36|l replicate the chunk 
medians to with 4%. Relative to the medians, while the 
optical chunks are relatively coherent, one can see the 
moving towards redder wavelengths causes the chunks to 
vary more about the median, for the reasons outlined 
above. 

3.3. Surface Gravity: BT-Settl Based 

Surface gravity, through pressure broadening effects, 
can change both the depth and width of spectral lines, 
and thus the RV velocity uncertainty of a spectrum. 
The exact response of lines to changes in surface grav¬ 
ity is dependent upon several factors. For example, a 
decrease in gravity can cause a weak line to either gain 
or lose strength depending upon the ionization state of 
the atoms. We therefore approached surface gravity ef¬ 
fects in a manner similar to effective temperature, by 
determining a purely numerical and descriptive scaling 
of velocity uncertainty. 

We took BT-Settl spectra from 2600K to 7600K with 
surface gravities of log(g) = 4.0 and log(( 7 ) = 5.0 and 
calculated the velocity uncertainties of each relative to 
log( 5 ) = 4.5. (Figure |T3|). We again split each spec¬ 
trum up into 100 A chunks and examined all three of 
our wavelength bands. Figure [13] displays the results for 
the individual chunks, as well the median velocity uncer- 



Fig. 14.— Th e points show the measured median slopes used 
in Equation m to describe the change in velocity precision as a 
linear function of the change in surface gravity (logg). The over- 
plated lines are the fits to these measurements, as per Equation 


tainty across all the chunks, and our fit to the median 
as a function of effective temperature and surface grav¬ 
ity. All three bands behaved roughly the same, with the 
greatest relative difference occurring around 4500K and 
the smallest differences happening towards the ends of 
the temperature range. 

Since the chunk medians in all three bands are nearly 
symmetric about unity, we decided to describe the 
change in velocity error as a linear function of surface 
gravity relative to what the velocity uncertainty would 
be for logg = 4.5. Specifically, 

/(log g) = m- A log(g) + 1 (37) 

with A log( 5 ) = logg — 4.5. The slope m depends upon 
the effective temperature and the band observed, 

TOOpt = -0.27505 (1 - 1.22211 ATe// (38) 

-4.17622 AT^^y) 

mRed =-0.33507(1 - 1.41362 AT^jf 

-4.63727 AT^^) 

=-0.43926(1- 1.12505 AT^ff 

-4.53938 AT^^), 

where ATg// = Te///5800K— 1. Figure [Ml shows the 
measured values of m across the temperature range, over¬ 
plotted by the fits from Equation (1551) . Note that these 
results are only accurate for surface gravities between 
4.0 < logg < 5.0. 

Since the extremes of our considered temperature 
range, where we expect |A log(( 7 )| to be the largest, show 
the smallest effect due to changing gravity, we expect 
surface gravity to play a relatively minor role in deter¬ 
mining the velocity error of a spectrum. Indeed, when 
we include how we expect surface gravity to vary with 
stellar mass (Equation [l^) we find that changing surface 
gravity is never more than a 10 % effect, and is frequently 
less. 

3.4. Metallicity: Kurucz92 Based 

Finally, we consider the effect of differing metallicity, 
which we parameterized as [Fe/H], on velocity uncertain¬ 
ties. Similarly to the effect of changing effective tem¬ 
perature and changing surface gravity, we expected that 
changing metallicity would alter the velocity uncertainty 
in a spectrum in a way that is difficult to capture ana¬ 
lytically from a priori arguments. We therefore confined 
ourselves to a numerical, descriptive, scaling for the effect 
of changing [Fe/Hj. 

To do so, we used three sets of Kurucz92-model spec¬ 
tra at fixed effective temperatures of 5000K, 5750K, 
and 6500K. All three sets had a fixed surface grav¬ 
ity of logg = 4.5, no a-enhancement, and metallicities 
of [Fe/H] = [-2.0,-1.5,-1.0,-0.5,+0.0,+0.5]. We used Ku- 
rucz92 based spectra, instead of BT-Settl spectra, due 
to the availability of Kurucz92 models with a wide range 
of [Fe/H] values. The BT-Settl spectra, while more phys¬ 
ically motivated in terms of what metallicities are avail¬ 
able for a given effective temperature and surface gravity, 
do not provide an arbitrary range of metallicities with no 
a-enhancement. As a result, our results are confined to 
the optical wavelengths, from 4000 A to 6500 A, where 
the Kurucz92 line-lists are robust. 




Beatty & Gaudi 


14 



[Fe/H] 


Fig. 15.— The velocity uncertainty of 100 A chunks as a func¬ 
tion of metallicity using Kurucz92-based spectra. These results are 
for a star with a fixed effective temperature of 5750K and a fixed 
log g = 4.5 in the optical (4000 A to 6500 A). The top panels shows 
the absolute uncertainty for the chunks, while the middle row shows 
the velocity uncertainty for each chunk normalized to [Fe/H]=0.0. 
The red points in the middle panel are the median of the relative 
velocity uncertainties. The bottom row shows the fractional dif¬ 
ference between each individual chunk’s relative uncertainties and 
this median. As expected, a decrease in [Fe/H] reduces the number 
of lines and average line strength in a spectrum, thus increasing the 
velocity uncertainty. 


Figure [15] shows the results of varying the metallicity 
on the velocity uncertainties for a 5750K, logg = 4.5 
star in the optical. The results for the 5000K and 6500K 
spectra were similar, differing by at most 15% at the 
low metallicity end. As one would expect, the velocity 
uncertainties for all of the chunks increases as [Fe/H] de¬ 
creases; a result of the absorption features in the spectra 
becoming weaker and less numerous. Some of the spe¬ 
cific 100 A chunks are strongly effected by this, as their 
major lines are drastically diminished at [Fe/H]=-2.0. 

We fit to the red chunk median in the middle panel 
of Figure [TS] to find the change in velocity uncertainty 
relative to [Fe/H]=0.0. This gave 

/([Fe/H]) = (39) 

For all three effective temperatures, this result is accurate 
to 15% over the range of [Fe/H] values we considered, 
with the highest difference occurring at [Fe/H]=-2.0. 


3.5. Final Expressions 

Putting together all of our results from the preceding 
analysis, we arrive at a semi-analytic expression for the 
velocity uncertainty using an arbitrary number of 100 A 
chunks as a function of R, usini, T^ff, logg, and [Fe/H]: 


1 
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Fig. 16.— The radial velocity uncertainty (in km s“^) for indi¬ 
vidual 100 A chunks in the optical, over the temperature range we 
consider. This assumes a conti nuu m flux in each chunk equal to 
unity, or in terms of Equation (I40II . that Jo,i = 1 photons per km 
s“^. This is a display of the information in Table 1. The darker 
regions indicate a lower velocity uncertainty at a given wavelength 
and effective temperature. 



VO.21660^ + 0^ +0.518 + 

00 (Te//) 

xfiT^ff) /(logs) /([Fe/H]). 

Where 00(7///) is given in Equation (|3T|) . f{Teff) in 
Equation (13^ . /(logs) in Equation (l37l) and /([Fe/H]) 
in Equation (|39|1 . The leading summation term is a sum 
over all of the 100 A chunks observed, with ^ as the 
velocity uncertainty of the individual chunks for R=oo, 
Vrot=Q, Vmac=0, and for a continuum level of 1 photon 
per velocity element. The true continuum level in each 
chunk is incorporated via Note that we are defin¬ 

ing the continuum here in velocity-space, and not pixel- 
space as is conventional in the observational literature. 

Table 1 lists values of in km s“^, normalized to 
Io,i = 1 photon per velocity element, for 100 A chunks 
between 4000 A and 25000 A for temperatures between 
2600K and 7600K and a constant logs=4.5. This infor¬ 
mation is also displayed graphically for our three wave¬ 
length regions, in Figures [TB] [T3 and [TB] Table 1 was 
calculated using the BT-Settl spectra, which used the 
lAsplund et al.l (1200911 solar abundances and had solar 
metallicity with no a-enhancement. The values in Ta¬ 
ble 1 are normalized for a continuum level of 1 photon 
per km s“^. Spectroscopic observations usually quote 
their SNR per pixel {SNRpi^^i), which can be converted 
into the appropriate units for /oy by taking 


Io,i — 


{SNRpi^^i)'^ Up 


AW 


chunk.i 


(41) 


where AVchunk,i is the velocity span of the wavelength 
chunk, and npix,i is the number of pixels on the detector 
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TABLE 1 

Chunk Velocity Uncertainties (km/s) 



2600K 

2800K 

3000K 

3200K 

3400K 

3600K 

3800K 

4000K 

4200K 

4000 A 

0.0681 

0.0756 

0.0821 

0.0919 

0.0930 

0.0856 

0.0911 

0.0821 

0.0758 

4100 A 

0.0783 

0.0829 

0.0848 

0.0916 

0.0873 

0.0821 

0.0739 

0.0655 

0.0556 

4200 A 

0.1232 

0.1175 

0.1236 

0.1413 

0.1387 

0.1302 

0.1186 

0.1059 

0.0887 

4300 A 

0.0487 

0.0597 

0.0769 

0.0930 

0.1062 

0.1106 

0.1109 

0.0999 

0.0923 


Note. — The complete table is available as online data at http://www.personal.psu.edu/tgbl5/beattygaudi/tablel.dat 
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stars cooler than 7600K: 


Te/f = 5603K ( ^ 
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Fig. 17.— The radial velocity uncertainty (in km s~^) for in¬ 
dividual 100 A chunks in the red, over the temperature range we 
consider. This assumes a conti nuu m flux in each chunk equal to 
unity, or in terms of Equation Hot . that /o,i = 1 photons per km 
s“^. This is a display of the information in Table 1. The darker 
regions indicate a lower velocity uncertainty at a given wavelength 
and effective temperature. 

used to observe the chunk. 

4. STELLAR PROPERTIES 

Having so far considered the dependence of velocity un¬ 
certainties on stellar parameters independently of what 
is physically reasonable, we know wish to apply this for¬ 
malism towards real stars. Our focus here is to determine 
what is the dominant source of the velocity uncertainty 
for main-sequence stars. 

We have chosen to focus on main sequence stars so 
that we may use stellar mass as a single variable to then 
calculate all of the stellar properties that determine ve¬ 
locity uncertainty. As above, the stellar properties we 
are interested in are rotation velocity, macroturbulence, 
effective temperature, and surface gravity. Additionally, 
we also wish to know the overall bolometric luminosi¬ 
ties and radii of the stars. The first allows us to calcu¬ 
late the continuum level of the spectra, while the second 
will be necessary to estimate rotation velocities from the 
rotation periods predicted from stellar gyrochronology 
relations. To that end, we fit relations for effective tem¬ 
perature, luminosity, radius, and surface gravity from the 
measurements listed in Table 1 of iTorres et al.l (|201Clll for 



For all the functions except the radius relation we have 
also given rough approximations, which we use when sim¬ 
plifying our numeric results. We note that these relations 
are only roughly consistent with each other, a result of 
our collapsing stars of different main sequence ages and 
metallicities onto single relations. Unlike temperature, 
luminosity, radius, and surface gravity, the rotation ve¬ 
locity does not display a simple scaling with stellar mass, 
and so we consider it separately and in more detail in the 
next section. 

For the macroturbulent velocity, we used the empiri¬ 
cal relation ior Vm.ar function of effective temperature 
determined bv iBruntt et al.l (|2010D . 

= 1.976 km s“^ -b 16.14 ATe// -b 19.713 

(46) 

where we have sub stituted ATgff = Te///5800K — 1. 
iBruntt et al.l (|2010ll make the point that this relation is 
only valid for stars with log((7) > 4.0 and between 5000K 
to 6500K. We use Equation (l46ll to estimate the macro- 
turbulent velocity for stars from 5000K to 7600K, and set 
a constant macroturbulent velocity of 0.51km (the 
value of Equation (H5)) at 5000K) for all stars cooler than 
5000K. Above 6500K, thus should not introduce large er¬ 
rors into our results, because the rotational velocities of 
these hot stars are at least five times larger than the 
calculated macroturbulent velocities. 

In addition to these stellar properties, there is also RV 
“jitter” in stars, which causes additional uncertainty in 
precision velocity measurements. Astrophysical jitter is 
generally a result of either star-spots on the stellar pho¬ 
tosphere, or short-period solar-like asteroseismic oscil¬ 
lations. For the latter, many of the existing RV sur¬ 
veys mitigated the effect of asteroseismic jitter by in¬ 
tegrating on a star for longer than the oscillation peri¬ 
ods, which are typically about 5 minutes. Unfortunately, 
the jitter caused by star-spots has no comparable solu¬ 
tion, other than avoiding stars with high activity indices. 
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Fig. 18.— The radial velocity uncertainty (in km s~^) for in¬ 
dividual 100 A chunks in the NIR, over the temperature range we 
consider. This assumes a conti nuu m flux in each chunk equal to 
unity, or in terms of Equation llioll . that 7o,i = 1 photons per km 
s~^. This is a display of the information in Table 1. The darker 
regions indicate a lower velocity uncertainty at a given wavelength 
and effective temperature. The sharp feature around 23000 A is 
the CO molecular bandhead. 

Typical jitter values for main-sequence stars cooler than 
6300K and with averag e activity levels are measured to 
be around 3 to 4 m s~^ (IWrightll2005Hlsaacson fc Fischer! 
l201Ct iMartmez-Arnaiz et alJ 12010(1 . As has been noted 
by all of these authors, this undoubtedly includes jit- 
ter from astrophysi c al and instrumental sources. Indeed, 
llsaacson fc Fischeil (|2010fl find that their measured jitter 
in K-dwarfs is completely uncorrelated with stellar activ¬ 
ity; they therefore conclude that the jitter displayed by 
these stars is likely a result of instrumental effects. 

Since velocity jitter is caused by additional astrophysi¬ 
cal and instrumental sources, the appropriate way to in¬ 
corporate it into our formalism is to add the jitter value 
for a given star in quadrature to the Poisson velocity 
uncertainty calculated using Equation (HOI) . RV jitter is 
therefore not a component we need to consider in deter¬ 
mining the dominant sources of Poisson velocity uncer¬ 
tainty in main-sequence stars. We thus leave it aside for 
now, other than to note the importance of jitter in using 
our results to fully model a realistic RV survey. 

4.1. Stellar rotation 

Since stellar rotation strongly broadens stellar lines, we 
undertook a detailed examination of the true rotational 
speeds, Vrot, of stars within our mass range. In general, 
stars with substantial outer convective envelopes, from 
O. 4 M 0 to the Kraft Break (|KraftJ FlQTflf l at 1.1 Mq will 
magnetically brake over the first billion years of their 
lives and coalesce onto a single mass-rotation-a ge rela¬ 
tion. This is the basis of stellar gyrochronology (iBarnesI 
1200311 . Stars less massive than 0.4Mq generally do not 
brake effectively, and so do not evolve onto a single mass- 
rotation-age relationship. Similarly, stars more massive 


than the Kraft Break mass of 1.1 Mq have very thin outer 
convective envelope and retain almost all of their pri¬ 
mordial angular momentum. These stars will slightly 
lengthen their rotational periods due to the gradual in¬ 
crease of their radii on the main sequence, but this is 
change is on the order of 2% over their lifetime. These 
heavier stars also do not, therefore, evolve onto a single 
mass-rotation-age relationship. 

For stars between 0.4Mq and I.IMq we used the 
modified Kawaler spin dow n model developed by 
Epstein & Pinsonneaultl (j2014f l to determine the rotation 
periods of stars at a certain mass and age. We then used 
the mass-radius relation in Equation (I44p to convert the 
rotation periods into rotation speeds. For stars older 
than 0.5 Gyr this spin-down model predicts a tight mass- 
rotation-age relation down to O. 4 M 0 , with more scatter 
as one goes to younger ages and lower mass. We linearly 
interpolated between the available model grid points in 
mass and age and took the median rotation period as the 
rotational period of all the stars with that mass and age. 

For stars less massive than O. 4 M 0 , we treated Vrot as a 
distribution, with velocities uniformly distributed in ve¬ 
locity between zero and some upper bound z^max- This 
roughly replicates the distribution of M-d warf rotation 
velocities observed bv iReiners et al.l (|2012ll . The upper 
bound was set equal to the Kawaler rotation velocity at 
the high mass end, Vrot(0.4MQ), which is an age depen¬ 
dent quantity, and increased linearly with mass through 
10 km s“^ at 0.2 Mq. Thus 


Vvot(M^) = ( [ 10 km s ^ - VrotiHAMo)] 


O.4M0 - M, 


(47) 


0.2Mq 


+ Vrot{0-‘1MQ) 

for M* < O. 4 M 0 . 


We treated stars more massive than I.IM 0 in a simi- 
lar manner. Base d on the observations and discussion in 
iGaigd (|1993ll and IReiners fc SchmitO (|2003ll , we treated 
the V sin i distribution of stars heavier than 1. IM 0 as uni¬ 
formly distributed in velocity between zero and a mass 
dependent upper bound. For the massive stars, this up¬ 
per bound was set to the Kawaler rotation velocity at 
I.IM 0 , Vrot{l.lMQ), and the bound increased linearly 
with mass through 100 km s“^ at 1 . 5 M 0 . Therefore for 
massive stars we have 


(M,)= [100km s-^-z;^ot(l.lM 0 )] (48) 


M* - I.IM0 


0.4Mq 


- Vrot{l.lMQ) 

for M, > 1 . 1 M 0 . 


5. DOMINANT SOURGES OF VELOGITY 
UNGERTAINTY 

We now wish to assess the dominant astrophysical 
sources of velocity uncertainty in RV measurements of 
main sequence field stars. To do so, we calculated the 
change in the velocity uncertainty caused by the /o,i, 0 , 
Teff, logg, and [Fe/H] terms in Equation (|i(Il) . This 
gives us the relative change in velocity uncertainty as a 
function of stellar properties, independent of the specific 
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wavelength chunk (or chunks) chosen. The specific wave¬ 
length information is provided in Equation (1401) by the 
ay values given in Table 1, and serves to simply set the 
appropriate absolute value of the uncertainty. 

We begin by calculating how the relative velocity un¬ 
certainty scales as a function of stellar mass in the wave¬ 
length ranges we consider. We used the relations in Sec¬ 
tion 4 to determine the effective temperature, surface 
gravity, and luminosity as a function of mass. For the 
luminosity, we included the effect of overall changes in 
the bolometric luminosity, normalized to a 5800K star, 
by Equation (1^ . and the effect on the observed lumi¬ 
nosity caused by the shifting of the blackbody emission 
across the specific wavelength range being observed. We 
refer to this as the “blackbody effect.” This blackbody 
term means that the exact results will still depend on the 
specific wavelength chunk used to calculate the relative 
uncertainties. For each of our three bands (optical, red, 
and NIR), we used a 100 A chunk in the middle of the 
wavelength range to calculate the blackbody effect. In 
our tests, using chunks at the extreme of our wavelength 
bands changes the calculated uncertainties by 5% or less, 
and does not affect our ultimate conclusions. 

As a fiducial example, we set i?=60,000, [Fe/H]=0.0, 
and the stellar age to 2.0 Gyr. Recall that the age 
will set the rotation velocity of stars between 0.4 Mq 
and 1.1 Mq, with younger stars rotating more rapidly. 
We chose 2.0 Gyr so as to be broadly representative of 
a typical field FGK d warf in the Solar neighborhood 
(|Nordstr6m et al.l 1200411 . 

Figures [IHl Uni and \Tn show the expected uncertainty 
as a function of mass for the optical, red, and NIR as the 
solid red line. To illustrate how the overall uncertainty 
is determined by the underlying stellar parameters, these 
figures also show how the uncertainty changes if we fix 
all but one of the physical processes that affect the ve¬ 
locity uncertainty and depend on stellar mass. These 
parameters are changes to line strengths and numbers 
due to effective temperature (T), changes in the overall 
bolometric luminosity (L), luminosity changes from the 
blackbody peak shifting relative to the spectral bandpass 
(BB), stellar rotation (VS), macroturbulence (VM), sur¬ 
face gravity changes (G), and spectral resolution (i?). 

There are two things to immediately note. First, below 
O.4M0 and above I.IMq the VS curve is for the maxi¬ 
mum observed usini at each mass. Second, the changing 
effect of spectral resolution as a function of mass is a re¬ 
sult of the average inherent line widths varying with ef¬ 
fective temperature, per Equation m- Larger inherent 
line widths (e.g., at lower temperatures in the optical) 
cause finite spectral resolution to have a smaller effect 
on the velocity uncertainties. 

In general, these three figures demonstrate that there 
are two general regimes for the RV errors of F-M main 
sequence stars: luminosity and temperature dominated 
uncertainties below I.IMq when stellar rotation is low, 
and rotation dominated uncertainties for stars above that 
mass. In particular, the rapid increase in the average 
rotation for more massive stars means that in all three 
bands the maximum velocity uncertainty raises sharply 
in this regime, becoming an order of magnitude larger 
than it would be for a Sun-like star at just ~ I.25Mq. 

Below I.IMq, on the other hand, the three wavelength 



Fig. 19. — The velocity uncertainty of observations of “realistic” 
stars, using the stellar properties described in Section 4, as a func¬ 
tion of stellar mass in the optical. The uncertainties are shown rela¬ 
tive to the uncertainty of an observation of a “perfect” non-rotating 
Sun-like star with no macroturbulence and observed using i?~oo. 
The red line shows the exact uncertainty calculation for observa¬ 
tions using a spectral resolution of /?=60,000 and a stellar age of 
t=2.0 Gyr. We have decomposed the overall velocity uncertainty 
into its constituent parts (black lines), as described in the third 
paragraph of Section 5. These are: changes to line strengths and 
numbers due to effective temperature (T), changes in the overall 
bolometric luminosity (L), luminosity changes from the blackbody 
peak shifting relative to the spectral bandpass (BB), stellar ro¬ 
tation (VS), macroturbulence (VM), surface gravity changes (G), 
and spectr al r esolution (R). The blue approximation line is from 
Equation Hot . and has been normalized to match the calculated 
uncertainty at 0.5Mq. In the optical, the approximated uncer¬ 
tainty is roughly proportional to the overall luminosity (L). 

regions behave differently. This primarily due to the 
changing effect of the blackbody peak shifting relative 
to the observed bandpass (the BB line in all three fig¬ 
ures). In the optical (Figure [121) the blackbody effect 
causes larger uncertainties for lower mass stars as the 
peak shifts into the red, particularly below O.GMq. In the 
red (Figure [20]), the blackbody term is nearly constant, 
as the peak is moving through this wavelength regime. 
By the time we reach the NIR fFigure I^Tl) . the blackbody 
term finally begins to reduce the velocity uncertainties of 
lower mass stars relative to solar-mass stars. 

To illustrate this different behavior in the different 
bandpasses, and to understand the dominant source of 
velocity uncertainty in a conceptually straightforward 
manner, we derived simple approximations to how the 
velocity uncertainty of low-mass (M* < I.IMq) stars 
scaled with stellar mass. In the optical, the blackbody 
effect nearly cancels the effect of changing effective tem¬ 
perature, leaving the bolometric luminosity as the dom¬ 
inant source of uncertainty. In the red, where the black¬ 
body effect is nearly constant, we approximate the ve¬ 
locity uncertainty as a combination of luminosity and 
temperature effects. Finally, in the NIR, we must also 
include the blackbody effect into our approximation us- 
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Fig. 20. — The velocity uncertainty of observations of “realis¬ 
tic” stars, using the stellar properties described in Section 4, as a 
function of stellar mass in the red. The uncertainties are shown 
relative to the uncertainty of an observation of a “perfect” non¬ 
rotating Sun-like star with no macroturbulence and observed us¬ 
ing Rr^oo. The red line shows the exact uncertainty calculation for 
observations using a spectral resolution of /?=60,000 and a stellar 
age of r=2.0 Gyr. We have decomposed the overall velocity un¬ 
certainty into its constituent parts (black lines), as described in 
the third paragraph of Section 5 and in the ca ptio n of Figure 19. 
The blue approximation line is from Equation (SSJ, and has been 
normalized to match the calculated uncertainty at 0.5 Mq. In the 
red, the approximated uncertainty is roughly proportional to the 
overall luminosity (L) and temperature effects (T). 

ing luminosity and temperature. We find that in the NIR 
the BB term is well fit by a simple linear relation with 
stellar mass, such that (Tv,bb oc OA^M^/Mq) + 0 . 6 . 

For our approximations for the relative uncertainty 
scaling below I.IMq we therefore have: 

/L 

cry,Opt oc ( ^ j for M* < I.IMq (49) 

/L 

cry, Red oc fRed{Teff) f j for M* < I.IMq 

/L 

Cry,NIR oc [O. 4 {M*/M 0 ) + 0.6] fNIR{Teff ) 

for M* < I.IMq. 

Above I.IMq, where stellar rotation dominates, we 
would approximate the velocity error as cry oc oo, as the 
rapidly increasing rotation velocities widen the lines to 
unusability for detecting all but the most massive com¬ 
panions. Note that, in detail, since values of usini are 
observed to be uniformly distributed in this high-mass 
regime, some stars will have a v sin i low enough that 
rotation will not dominate the velocity uncertainties. 

The red lines in Figures 15, 16, and 17 show these 
approximations. Note that we have normalized each 
to match the directly calculated velocity uncertainty at 
0.5Mq. This serves to account for the added uncer¬ 



Fig. 21.— The velocity uncertainty of observations of “realis¬ 
tic” stars, using the stellar properties described in Section 4, as a 
function of stellar mass in the NIR. The uncertainties are shown 
relative to the uncertainty of an observation of a “perfect” non¬ 
rotating Sun-like star with no macroturbulence and observed us¬ 
ing Rr^oo. The red line shows the exact uncertainty calculation for 
observations using a spectral resolution of R=60,000 and a stellar 
age of r=2.0 Gyr. We have decomposed this overall velocity un¬ 
certainty into its constituent parts (black lines), as described in 
the third paragraph of Section 5 and in the ca ptio n of Figure 19. 
The blue approximation line is from Equation ll49ll . and has been 
normalized to match the calculated uncertainty at 0.5Mq. In the 
NIR, the approximated uncertainty is roughly proportional to the 
overall luminosity (L), temperature effects (T), and the blackbody 
effect (BB). 

tainty arising from our assumed value for spectral res¬ 
olution, i?, and the non-zero rotation of the low-mass 
stars. One can see the approximate velocity uncertainty 
agrees well with the calculated velocity uncertainty for 
lower mass stars, but the two diverge as mass increases 
above M* rs O.SMq. This is a result of our ignoring the 
effects of stellar rotation in making our approximations. 

6 . DISCUSSION 

In addition to the above approximations for the scal¬ 
ing of velocity uncertainty as a function of stellar mass, 
and the general scaling of velocity uncertainty with stel¬ 
lar parameters given in Section 3, there are three general 
points regarding RV surveys that are interesting to con¬ 
sider. 

First, velocity precision does not scale linearly with 
vsini for high rotation, or as R~^ for low reso¬ 
lution, as has been claimed in s ome of the litera¬ 
ture (e.g.. iHatzes fc Cochranl Il99li iCoiines et al.l[T996l : 
IBouchv et al.l 120011: [Bottom et al.l I20I3I) . Instead, as 
we have shown, it goes appropriately as (wsinz)^/^, or 
j^- 3 / 2 ^ This arises due to two competing effects. First, 
the intrinsic velocity information in the line scales with 
the FWHM 0 of the line as 0^. Second, for weak lines, 
the photon noise is dominated by the continuum emis¬ 
sion. The continuum spanned by the line is oc 0, and 
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thus the photon noise scales a 0 The net result of 
these two effects results in the velocity uncertainty scal¬ 
ing as 0^/^. This is in contrast to an isolated line with 
negligible continuum, where the precision scales simply 
as 0. Additionally, it is also important to remember that 
the velocity uncertainty is set by the overall line width, 
which is a combination of the inherent line width and 
the effects from broadening mechanisms like v sin i and 
R. By considering the total resulting line width, we are 
thus able to fit for and describe the velocity precision as 
a function of v sin i and R over the entire range of these 
two parameters, and not just at the high usini or low R 
limits. 

Second, Figure H] illustrates how arbitrarily increasing 
the spectral resolution if an instrument does not give 
arbitrarily low velocity uncertainties. Specifically, there 
is a “knee” in the uncertainty curves around i?=60,000, 
after which increased spectral resolution has a much di¬ 
minished effect. This occurs because once a spectrograph 
reaches i?=60,000, the instrument is able to resolve al¬ 
most all the lines in a stellar spectrum; further increases 
to R therefore do not provide substantially more infor¬ 
mation. The precise spectral resolution we estimate for 
th e position of the “kn ee” is similar to that determined 
by iBouchv et al.l (j2001ll ,but h igher than the point of di¬ 
minishing returns estimated in iBottom et al.l (1201311 , who 
give a position of i?=45,000. This difference is a result of 
the relatively co arse wavelengt h samp ling in the BT-Settl 
spectra used by [Bottom et al.l (1201311 , which imposes an 
effective resolution floor of i?Ril 00 , 000 , as we discuss in 
more detail in Section 3 and around Figure 5. 

For stellar RV surveys. Figures [El [20l and |2T] also show 
that there is a limited utility to spectral resolutions above 
i?=60,000. One can see in these three figures that the 
effect of having R=60,000 (the horizontal dashed-line la¬ 
beled “R”) is the dominant source of velocity uncertainty 
for an extremely small range of masses. 

Instead, the dominant error source in radial veloc¬ 
ity measurements will either be caused by luminosity 
(M* < O. 8 M 0 ) or rotational velocity {O.SMq < M*). 
This means that relatively small detector arrays can be 
used effectively for multi-object RV surveys: 1024 pix¬ 
els along the spectral dispersion axis would allow for a 
100 A spectral chunk to be imaged at i?=60,000. As an 
example, if we ignore systematic and instrumental uncer¬ 
tainties, if such a survey observed a Sun-like star from 
5100-5200 A with SNR =200 per pixel, we predict that 
the Poisson velocity uncertainty on an individual obser¬ 
vation would be about 6 m s“^. 

More generally, we note that in terms of statistical ve¬ 
locity uncertainty there are sharply diminishing returns 
to be made from arbitrarily increasing spectral resolu¬ 
tion. First, above a certain resolution all of the ma¬ 
jor lines in a spectrum will be resolved; for resolutions 
beyond this value the dominant source of uncertainty 
will be set by the line widths themselves. From Equa¬ 
tion dSI]), we can see that this will happen for a Sun¬ 
like star for i?Ri60,000 in the optical, and i?Ri80,000 
and i?Ri45,000 in the red and NIR, respectively. Sec¬ 
ond, stellar rotation will limit the utility of increased 
spectral resolution - even in slowly rotating field stars. 
For i?=60,000, for example, the effective line broaden¬ 
ing caused by limited spectral resolution is equal to the 


amount of line broadening caused by 4.25 km s“^ of stel¬ 
lar rotation. For stars rotating faster than this, further 
increasing R will thus provide a small change in the mea¬ 
sured velocity uncertainty. Generically, the limiting spec¬ 
tral resolution for a given stellar rotation velocity will be 

Riim _ 4.25 kms“^ 

60,000 “ Vrot ■ 

That being said, extremely high spectral resolutions may 
help with correcting for systematics in observations, by 
enabling one to study the detailed shapes of the lines. 

Third, Figures [El (13 and [T51 show the importance of 
choosing the appropriate wavelength for an RV survey, 
particularly if one cannot cover a wide range of wave¬ 
lengths. To investigate this in more detail, we collapsed 
the uncertainty values behind Figures [El [El [El and Ta¬ 
ble 1 along the temperature axis to see what are the best 
locations for observations. We divided the temperature 
range into M-stars (2600K to 4000K), K-stars (4000K to 
5200K), G-stars (5200K to 6000K), and F-stars (6000K 
to 7600K) and took the median uncertainty values across 
these ranges for each 100 A wavelength chunk in Table 1 . 
Figure [22l shows the results, visualized in four different 
ways. First, the upper left panel shows the median uncer¬ 
tainties of each 100 A chunk for the four spectral types 
(labeled “Raw (linear)”). The lower left panel shows 
these same chunk medians, but now we have normalized 
them according to the fraction of the overall blackbody 
luminosity that each chunk occupies (labeled “BB Nor¬ 
malized (linear)”). In terms of Figures [El [SHI and [U 
this factors in the “BB” line. The two right panels in 
Figure [22] are similarly “Raw” or “BB Normalized,” but 
we have now combined the chunks into equally wide loga¬ 
rithmic bins, rather than plotting them linearly as before. 
This is meant to replicate the true observing mode of 
spectrographs: at fixed resolution an instrument can im¬ 
age proportionally more of a spectrum at proportionally 
longer wavelengths. The lower right panel of Figure [22] 
therefore most directly informs the selection of a proper 
observing wavelength. 

In the optical, spectra longwards of 5500 A provide rel¬ 
atively little information for F-stars, and provide an in¬ 
ferior amount of velocity information for later spectral 
types as compared to most of the shorter wavelengths. 
For a survey covering all spectral classes with a limited 
wavelength range, the best wavelengths to look at in the 
optical would be from 5000 A to 5200 A, which contains 
the M g b triplet, and has been noted previously (iLathaml 
Il985f) . Shorter wavelengths than this, while advanta¬ 
geous for FG stars, provide little gain for either K- or 
M-stars. 

In the “red,” between 6500 A and 10000 A, the best 
regions for a general exoplanet survey are in the I or i' 
bands around 7000 A. Moving further out into the red, to 
z', provides less velocity information at almost all effec¬ 
tive temperatures. Indeed, as shown in the bottom right 
panel of Figure [ 22 ] observing from 7000 A to 8000 A in I 
or i' has the lowest base uncertainty for observations of 
M-dwarfs over the entire wavelength range. 

To investigate this in more detail. Figure [E] shows four 
of the temperatures that make up the M-dwarfs in Fig¬ 
ure [E] displayed in a similar manner. One can see that 
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Fig. 22.— The base uncertainty values from Table 1 divided 
into M-stars (2600K to 4000K), K-stars (4000K to 5200K), G- 
stars (5200K to 6000K), and F-stars (6000K to 7600K) and me¬ 
dian combined within those categories. The two left panels show 
the medians for each 100 A wavelength chunk in Table 1, while the 
right panels show the median values after the chunks have been 
organized in 20 logarithmically-spaced bins. This is meant to repli¬ 
cate a spectrograph observing at fixed resolution. The top panels 
show the raw numbers from Table 1, while the bottom panels are 
normalized to according to the fraction of the overall bolometric 
luminosity that each chunk occupies. 

between 3800K (an MO) and 2600K (roughly an M7) the 
shape of the raw uncertainties as a function of wavelength 
are roughly self-similar, but the flux normalized uncer¬ 
tainties flatten out as one moves to cooler temperatures. 
This is, of course, a result of the peak of the stellar SED 
moving towards longer wavelengths as the temperature 
decreases. Nevertheless, although the peak of a 2600K 
blackbody has moved out to 11,000 A, one can clearly 
see in the lower panels of Figure [23] that the most ef- 
hcient wavelength region to observe all of the M-dwarf 
temperatures we consider is shortwards of 10,000 A, be¬ 
tween 6000 A and 9000 A. The M-dwarf spectra at these 
wavelengths have approximately 10 times the velocity 
information (i.e., 1/10 the uncertainty) as compared to 
the NIR, and thus remain the most efficient observing 
location even aft er accounting for the relative amount of 
stellar emission. iBottom et al.l (|2013[ 1 arrived at a simi¬ 
lar conclusion for the best wavelength regime to observe 
M-dwarfs. This illustrates the importance of considering 
the velocity information available in a wavelength region, 
and not just the shape of the stellar SED, when designing 
an RV survey targeted at M-dwarfs. 

If we now consider only the NIR, it is interesting to 
compare Figures [T51 and against the available trans¬ 
mission windows in the atmosphere. We can roughly 
approximate these windows by considering the NIR pho¬ 
tometry bands: J (11000 A to 13500 A), H (15000 A to 
17000 A), and K (20000 A to 23000 A). For a general RV 
survey of all spectral types, the optimum observing band 
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Fig. 23.— The base uncertainty values fr om Table 1 for four 
different M-dwarf temperatures. As in Figure[22l the two left pan¬ 
els show the uncertainties for the 100 A wavelength chunks listed 
in Table 1, while the right panels show the uncertainties after the 
chunks have been combined into 20 logarithmically-spaced bins. 
This is meant to replicate a spectrograph observing at fixed reso¬ 
lution. The top panels show the raw numbers from Table 1, while 
the bottom panels are normalized to according to the fraction of 
the overall bolometric luminosity that each chunk occupies. The 
bottom right panel particularly illustrates that the most efficient 
wavelength range to observe M-dwarfs is generally shortwards of 

10,000 A. 

would be in J, since the uncertainty for FGK spectral 
types steadily increases towards longer wavelengths. For 
a survey targeting only later spectral types, the decision 
is less clear cut. 77-band has the lowest uncertainty for 
observing M-dwarfs (down 10% compared to J and down 
40% compared to K). At the same time, if one wished to 
include K-dwarfs in this survey then the optimum observ¬ 
ing bandpass would again be J. It is interesting to note 
that the best wavelength regions to observe M-dwarfs 
are precisely in between J, H, and K. This is a result of 
molecular lines from water appearing in the cooler stars’ 
atmospheres, but this same water absorption in Earth’s 
atmosphere is precisely what sets the location of the NIR 
observing bands. 

7. SUMMARY 

We have considered the astrophysical sources of veloc¬ 
ity uncertainty in stellar RV measurements. In doing 
so, we are able to describe the basic mechanisms that 
cause velocity uncertainties, what the dominant driver 
is behind stellar velocity uncertainties at various stel¬ 
lar masses and in various wavelength regions, and fur¬ 
nish several points for consideration when designing an 
RV survey. In doing so, we emphasize that we have de¬ 
liberately decided to restrict our focus to astrophysical 
sources of uncertainty, and so we do not consider the 
effects of star spots, granulation, or asteroseismic pulsa¬ 
tions. Similarly, we do not treat instrumental velocity 
uncertainties like wavelength calibration, optical effects. 
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or instrumental drifts, nor do we consider the effect of 
telluric absorption lines on the spectra. 

We determine general scaling laws for the expected ve¬ 
locity uncertainty. This allows the reader to estimate the 
amount of velocity uncertainty present in observations of 
main-sequence stars using an arbitrary wavelength range 
between 4000 A and 25000 A, over a large set of possible 
spectral resolutions and stellar properties. This is in con¬ 
trast to previous work in this area, which has provided 
results using specific observing set-ups (e.g., fixed spec¬ 
tral resolution or stellar properties). 

At a basic level, we demonstrate that the velocity un¬ 
certainty of a weak spectral absorption line in a contin¬ 
uum scales as 0^/^, where 0 is the FWHM of the line, 
and not linearly with 0 as one expects when there is no 
continuum emission and as has been claimed in some pre¬ 
vious work. Using model spectra, we then calculated how 
the velocity uncertainty changes as a function of spectral 
resolution, stellar rotation, stellar effective temperature, 
stellar surface gravity, and stellar metallicity. By divid¬ 
ing our model spectra up into 100 A-wide chunks, we find 
that the effects of resolution, rotation, and surface grav¬ 
ity operate on the chunks in a largely self-similar manner 
- regardless of the specific wavelength or spectral fea¬ 
tures within a chunk. Effective temperature presents a 
more complicated picture, with different chunks behav¬ 
ing very differently. We numerically fit a rough relation 
to the chunk medians, but the variation between chunks 
as a function of temperature is one illustration of the im¬ 
portance of carefully choosing the wavelength range used 
in an RV survey. 

With these basic relations established, we are able to 
calculate how the velocity uncertainty scales as a function 
of stellar mass. For stars more massive than I.IMq, we 
find that the rapidly increasing stellar rotation dominates 
the predicted uncertainties. Below I.IMq, the velocity 
uncertainty is set by a combination of competing effects 
from changes in stellar luminosity, temperature and sur¬ 
face gravity. In the optical, between 4000 A and 6500 A, 
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we find that almost all of these effects cancel, leaving the 
velocity uncertainty to be predominately set by the bolo- 
metric luminosity of the target star for a fixed distance. 
This is not true in the red (6500 A to 10000 A) or the NIR 
(lOOOOA to 25000A), where one must also account for 
temperature (red) or temperature and the effect of the 
blackbody peak shifting relative to the observed wave¬ 
length range (NIR). We give simple approximations for 
how the velocity error scales with mass for each of these 
three wavelength regimes. 

More generally, our consideration of velocity errors in 
RV surveys highlights two important points for consid¬ 
eration. First, after a certain point increasing spectral 
resolution provides diminished returns. This primary oc¬ 
curs because once one has resolved the lines in a spec¬ 
trum increased resolution provides little more informa¬ 
tion, and because, depending on the stars being surveyed, 
stellar rotation will provide the dominant source of ve¬ 
locity uncertainty - not spectral resolution. We find that 
this point of diminishing returns occurs at approximately 
i?=80,000, though we note that extremely high spectral 
resolutions may help with correcting for systematics in 
observations. 

Second, the most efficient wavelength region to operate 
an RV survey for M-dwarfs is between 6000 A to 9000 A. 
Although the peak emission for M-dwarfs is generally 
longwards of these wavelengths, the base velocity uncer¬ 
tainties of spectra in this wavelength region are about 
1/10 that of spectra in the NIR bands. This means that 
even after accounting for the difference in received flux, 
M-dwarf spectra from 6000 A to 9000 A will give a lower 
velocity uncertainty than spectra in observed in the NIR 
at the same exposure time. 
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